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^ ■ cosmology and collider experiments on the allowed regions of parameter space as is 
'NT ■ 

I reduced from the GUT scale. We present separate analyses of the {mi/2,mQ) planes for 

^ ! tan/3 = 10 and tan/5 = 50, as well as a discussion of non-zero trilinear couplings, Aq. 

■ Specific scenarios where the gaugino and scalar masses appear to be universal below the 



GUT scale have been found in mirage-mediation models, which we also address here. We 
demand that the lightest neutralino be the LSP, and that the relic neutralino density not 



H I conflict with measurements by WMAP and other observations. At moderate values of M, 



we find that the allowed regions of the (mi/2,wio) plane are squeezed by the requirements 
of electroweak symmetry breaking and that the lightest neutralino be the LSP, whereas the 
constraint on the relic density is less severe. At very low M(„, the electroweak vacuum 
conditions become the dominant constraint, and a secondary source of astrophysical cold 
dark matter would be necessary to explain the measured relic density for nearly all values 
of the soft SUSY-breaking parameters and tan (3. We calculate the neutralino-nucleon cross 
sections for viable scenarios and compare them with the present and projected limits from 
direct dark matter searches. 
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1 Introduction 



Over the past three and a half decades, the Standard Model (SM) of particle physics has been 
remarkably successful at describing the interactions of elementary particles at or below the 
weak scale. However, there are several compeUing reasons to expect that the SM is merely a 
low-energy effective theory that fits into a larger framework. Chief among these reasons are 
the related hierarchy and naturalness problems, namely the creation and maintenance of a 
large hierarchy of mass scales despite the fact that the electroweak Higgs potential is unstable 
with respect to quantum corrections within the SM [1]. The appearance of supersymmetry 
(SUSY) at the TeV scale would not only solve the naturalness problem and facilitate the 
unification of gauge couplings at a high scale as in simple Grand Unified Theories (GUTs) [2], 
but also predict a light Higgs boson as apparently favoured by the high-precision electroweak 
data [3]. With the additional plausible assumption of R-parity conservation, the lightest 
SUSY particle (LSP) is stable and, if uncharged, is a natural candidate for astrophysical 
cold dark matter [4]. For these reasons, models with SUSY broken at the TeV scale are 
extensively studied. 

It is evident that SUSY must be broken, since we have not yet observed any superpartners 
of SM particles, but the mechanism of SUSY breaking and how this breaking is communicated 
to the observable sector have been the subjects of much speculation [5]. Phenomenologically, 
the magnitudes of the SUSY-breaking parameters observable at low energies are often cal- 
culated by assuming values of the soft SUSY-breaking parameters at some high input scale 
and evolving them down to lower scales using the renormalization-group equations (RGBs) 
of the effective low-energy theory. This is generally taken to be the minimal supersymmet- 
ric extension of the SM (MSSM) [6]. In the constrained MSSM (CMSSM) [7-13], the soft 
SUSY-breaking parameters are assumed to be universal at the high scale. It should be noted, 
however, that there are many theories of SUSY breaking in which the soft SUSY-breaking 
parameters are not universal at the input scale [14]. 

The CMSSM can be parametrized at the universality scale by five free input parameters, 
namely the scalar mass, mo, the gaugino mass, mi/2, the trilinear soft breaking parameter, 
Aq, the ratio of the Higgs vevs, tan/3, and the sign of the Higgs mass parameter, The 
input scale at which universality is assumed in CMSSM models is usually taken to be the 
SUSY GUT scale, Mqut ~ 2 x 10^^ GcV. However, it may be more appropriate in some 
models to assume the soft SUSY-breaking parameters to be universal at some different input 
scale. Mini which may either be intermediate between Mqut and the electroweak scale [15], 
the case studied here, or perhaps larger than Mqut [16]. 

Specific scenarios in which the soft SUSY-breaking parameters may be universal at a scale 
below Mgut occur in models with mixed modulus-anomaly mediated SUSY breaking (MM- 
AMSB), also called mirage-mediation [17], and models with warped extra dimensions [18]. 
In the case of mirage-mediation, the universality scale is the mirage messenger scale, which 
is predicted to be Mj„ ~ 10^° — 10^^ GeV in the case of KKLT-style moduh stabilization [19]. 
In other models, the universality scale may lie anywhere between 1 TeV and Mpi. 

In this paper, we present an in-depth study of the effect on the allowed regions of the 
CMSSM parameter space of lowering the assumed universality scale. We focus on the de- 
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pendences of the constraints from cosmology and collider experiments on the value of Mj„ 
in such GUT-less scenarios, paying particular attention to the regions of parameter space 
favored by the value of the cold dark matter relic density inferred from WMAP [20] and 
other measurements, assuming that the cold dark matter is mainly provided by the lightest 
neutrahno x- Within the GUT-less allowed regions, we also calculate the neutralino-nucleon 
cross sections and compare them with present and expected limits from direct searches for 
cold dark matter. 

This work is a sequel to the exploratory study of GUT-less CMSSM scenarios made 
in [15], in which our attention was restricted to the case tan/3 = 10, = 0, > and 
Min > 10^^-^ GeV. We found that, as the universality scale was reduced to this value, one of 
the most dramatic changes was to the footprint in the (mi/2, ttt-o) plane of the constraint on 
the relic abundance of neutrahnos inferred from WMAP et al. In the standard GUT-scale 
universality case, there are three well-defined cosmologically preferred regions of parameter 
space where the relic density of neutrahnos matches the estimate of the cold dark matter relic 
density based on data from WMAP and other observations: the coannihilation region [21], 
the rapid- annihilation funnel [9,22] and the focus-point region [23]. In the GUT-less GMSSM 
scenario [15], we found that, as the universality scale is lowered to Mj„ ~ 10^^ GeV, these 
regions approach and merge, forming a small WMAP-pref erred island in a sea of parameter 
space where the neutralino relic density is too small to provide all the cold dark matter 
wanted by WMAP. We found that, in this case, the only region with a neutrahno rehc 
density that exceeds the WMAP measurement is a 'vee' at large mi/2, bordering the region 
where the stau is the LSP. 

In this paper, we extend the previous analysis to include other values of Aq for tan j3 = 10, 
to the case tan/3 = 50, and to lower values of Mj„. For this purpose, we extend the code 
used previously to evaluate the cold dark matter density by implementing all coannihila- 
tions between the three lightest neutrahnos and the lighter chargino species. As we exhibit 
explicitly, their inclusion is essential for an accurate calculation of the relic density in some 
important regions of the GUT-less parameter space. The second objective of this paper is 
to calculate the neutralino dark matter scattering cross sections (both spin-dependent and 
spin-independent) in such GUT-less models. 

The outhne of the paper is as follows. In Section 2 we discuss briefly the renormal- 
izations of the SUSY-breaking contributions to the masses of the squarks, sleptons and 
gauginos as functions of Mj„, as a preliminary to provide background understanding for 
some of the results presented later. Then, in Section 3 we discuss the current experimental, 
phenomenological and cosmological constraints on CMSSM scenarios that we use. Section 
4 contains our core discussion of the variation in the allowed region of parameter space as 
Min is decreased from the GUT scale down to Mj„ = 10^ GeV, for both tan/? = 10 and 
tan/5 = 50. We also present a separate treatment of the mirage-mediation scenario [17] with 
KKLT moduli stabilization [19]. We then present in Section 5 the corresponding predic- 
tions for neutralino-nucleon scattering cross sections in GUT-less scenarios, and Section 6 
summarizes our conclusions. An Appendix motivates and discusses relevant details of our 
implementation of multi-channel neutralino and chargino coannihilation. 
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2 Renormalization of SUSY-Breaking Mass Parame- 
ters 



In order to understand the changes in the allowed regions in the {^1,1/2,1^0) plane of the 
CMSSM that occur as Mj„ is lowered, it is necessary first to understand the consequences 
for the observable sparticle masses of lowering the universality scale. In the CMSSM with 
universality imposed at the GUT scale, the one-loop renormalizations of the gaugino masses 
Ma, where a = 1,2,3, are the same as those for the corresponding gauge couplings, 
Thus, at the one-loop level the gaugino masses at any scale Q < Mqut can be expressed as 

MaiQ) = ""^if^ M MgutI (1) 

where Mai^McuT) = On the other hand, in a GUT-less CMSSM, where the gauge- 

coupling strengths run at all scales below the GUT scale but the soft SUSY-breaking pa- 
rameters run only below the lower universality scale, Mj„, at which all the gaugino masses 
are assumed to be equal to mi/2 = Ma{Min), we have 

at the one-loop level. Since the runnings of the coupling strengths in GUT and GUT-less 
CMSSM scenarios are identical, the low-energy effective soft gaugino masses, Ma{Q), in 
GUT-less cases are less separated and closer to mi/2 than in the usual GUT CMSSM, as 
seen explicitly in panel (a) of Fig. [T]0. 

The soft SUSY-breaking scalar masses are renormalized by both gauge and (particularly 
in the cases of third-generation sfermions) Yukawa interactions, so the running is somewhat 
more complicated. At the one-loop level one can summarize the effects of renormalizations 
at any Q < Mj„ as 

mliQ) = m2(M,„) + Q{Q,M,^)ml/^, (3) 

where we make the CMSSM assumption that the mQ(Mm) are universal at Mj„, and the 
Ci{Q, Min) are renormalization coefficients that vanish as Q ^ Min. We display in panels 
(b) and (c) of Fig. [1] the two-loop-renormalized soft SUSY-breaking masses of the the first- 
and second-generation left- and right-handed squarks, qL,R, the stop mass eigenstates, ti_2, 
and the left- and right-handed sleptons, El^r- We see again that in GUT-less cases the soft 
SUSY-breaking scalar masses are less separated and closer to mo than in the usual GUT-scale 
CMSSM. 

In the CMSSM, the electroweak vacuum conditions are used to fix the values of 
and mA- Although we use the full two-loop renormalizations, insight into the effects of 



^Note that in making this plot we have included the full two- loop renormalization- group equations for 
the gaugino masses, which are not identical to those for the gauge couplings, although the difference is not 
very striking. 
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Figure 1: The dependences of observable sparticle mass parameters on the input scale Mj„ at 
which they are assumed to be universal: (a) gaugino masses Mi^2,3; (b) squark masses, (c) 
slepton mases, and (d) Higgs (mi^2,mA), stau and the lightest neutralino x masses, as well 
as fi and the U(l) gaugino mass Mi. The calculations are made for the representative case 
mi/2 = 800 GeV, mo = 1000 GeV, Aq = 0, tan/5 = 10 and /i > 0. 
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varying M(„ on the required values of |/i| and rriA can be gleaned from simple leading-order 
expressions. The tree-level solution for /i is 



ml — ml tan^ (3 M| 
tan2 (3-1 2 



= ^ _ (4) 



where mi and m2 are the soft Higgs masses associated with Hi and H21 respectively. The 
variation of with M^^ for one fixed pair of values of (mi/2,'"io) is seen in panel (d) of 
Fig. [H where we see that the solution of (jl]) for /x^ becomes negative and unphysical for 
Min < 10^° GeV. For this value of M(„, values of mo > 1000 GeV would not yield physical 
electroweak vacua. One can see from ^ and panels (b) and (c) of Fig. [T] that, as Mj„ 
decreases, the soft scalar masses remain closer to the input value, mg. This has the converse 
result that, for any fixed mi/2, as the universality scale is lowered, /i^ changes sign and 
becomes unphysical at smaller values of mo, causing the upper boundary of the unphysical 
region to creep down farther into the (mi/2,'^0) plane. This explains the encroachment of 
the upper-left excluded regions in the (mi/2,mo) planes shown later in Figures [2] - El as M(„ 
decreases. 

The weak-scale value of m^ decreases with M^n logarithmically, as also seen in panel (d) 
of Fig. [H and also in panels (c) and (d) of Fig. 3 of [15]. In addition to its importance for 
the direct detection of the near-degenerate A, H and H"^ bosons, this feature is important 
indirectly for several aspects of our later discussion. One is the constraints from heavy- 
flavour physics to be discussed in the next section: since 6 — > 57 and Bg — > yU+ yU~ at large 
tan (3 have important contributions from the exchanges of heavier Higgs bosons, the impact 
of these constraints increases as m^ decreases and hence as Mj„ decreases. A second impact 
of m^ is on the cold dark matter density: since a rapid-annihilation funnel appears when 
m^ ~ for fixed values of the other parameters such as tan /5, mo and Aq, this funnel 

appears at lower m^ and hence mi/2 as Mj„ decreases. Finally, another potential impact 
is on the spin-independent neutralino dark-matter scattering cross section, which receives a 
significant contribution from heavy Higgs exchange, as discussed later. 

In addition to the excluded regions in the upper left corners of each of the (mi/2,mo) 
planes shown in Figures [210 where electroweak symmetry breaking is not obtained, we see 
a second major excluded region in the lower right corner of each panel. In these regions of 
the plane, the lightest stau, fi, becomes lighter than the lightest neutralino, resulting in a 
charged LSP, which is incompatible with general arguments from astrophysics and cosmology. 
As we see from ([3]), as decreases the positive coefficient Cf^ also decreases because Min 
is approaching the low scale, Q. Hence m^^ gets progressively closer to mo for any fixed mi/2, 
as seen in panel (c) of Fig. [H At the same time, the gaugino masses remain closer to mi/2 as 
Min decreases, implying that, as long as the lightest neutralino remains essentially a bino, its 
mass becomes a larger portion of the universal gaugino mass. This can be seen in panel (d) 
of Fig. [H where for this particular point in the (mi/2, '^o) plane, the LSP mass tracks that of 
the bino for Mj„ > 10^^ GeV. As a result, for fixed mi/2 and mo, as the universality scale Mj„ 
is lowered from Mquti initially m^^ increases and m^^ decreases. Hence, as Mj„ decreases 
for any fixed mi/2, a larger value of mo is required to enforce the condition m^^ < m^^. For 
this reason, the lower-right excluded regions in the (mi/2,mo) planes shown in Figures [2] - [5] 
initially expand as Min decreases. 
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However, since decreases as Mj„ decreases, as discussed above, below a certain value 
of Min, |/i| becomes small enough that the lighter Higgsino takes over as the LSP, with a 
mass that decreases as continues to decrease. In panel (d) of Fig. [1], one can see that, for 
Min ^ 10^^ GeV, the LSP is sufficiently Higgsino-like that its mass is nearly identical to 
Since the boundary of the disallowed stau LSP region is determined by equality between the 
masses of the stau and the lightest neutralino, this boundary therefore falls to lower rriQ when 
Min is decreased below the bino-Higgsino cross-over point, as is seen in in the (mi/2,'"io) 
planes shown later in Figures [2] - [51 

3 Experimental, Phenomenological and Cosmological 
Constraints 

Our treatments of experimental, phenomenological and cosmological constraints essentially 
follow those in [15], but with differences that we describe below. 

3.1 LEP Experimental Constraints 

The appropriate LEP lower limit on the chargino mass for the class of CMSSM models 
discussed here is m^± > 104 GeV [24], and the nominal effective lower limit on the mass 
of the lightest Higgs boson /i is 114 GeV @ [25]. However, in addition to displaying the 
direct position of the 114 GeV bound in the GUT- less parameter space, we also calculate 
and display the 95% CL limit obtained by combining the experimental likelihood, Cexp, from 
direct searches at LEP 2 and a global electroweak fit, convolved with the theoretical and 
parametric errors in rrih which provides a more exact (and relaxed) interpretation of the 
LEP Higgs limit within the MSSM. The top mass used in these calculations is rrit = 171.4±2.1 
GeV [26]. 

We note that one can use ([3]) to predict how the impact of the LEP Higgs mass constraint 
varies with Mj„. We recall that the mass of the lightest scalar MSSM Higgs boson ruh < Mz 
at tree level, but is renormalized by an amount that depends logarithmically on mj. Eq. ([3]) 
shows that mj- decreases as Min is lowered for any fixed mi/2 and itlq. However, mj- also 
increases with mi/2- Thus, one should expect that the LEP Higgs constraint moves to larger 
mi/2 as the universality scale is lowered. 

3.2 Muon Anomalous Magnetic Moment 

It is well known that the measurement by the BNL g-2 Collaboration [28] disagrees signifi- 
cantly with the Standard Model if e^e~ annihilation data are used to calculate the Standard 

^We implement this constraint by calculating the lightest Higgs mass with the previous version of the 
FEYNHIGGS code [27], which incorporates a direct interface with the underlying CMSSM parameters, and 
allowing a possible error of 1.5 GeV to account for possible higher-order contributions. We have verified that 
the numerical difference from the more recent version of FEYNHIGGS is considerably smaller than our error 
allowance. 

^ We thank A. Read for providing the LEP CLs values. 
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Model contribution, but there is no significant discrepancy if tliis is calculated using r-decay 
data [29]. In view of the lack of consensus on the interpretation of the measurement of 
O'fj. = {g^i — 2)/2, we use it only as part of our motivation for restricting our study to the case 
H > 0. However, if the e'^e~ estimate of the hadronic contribution to the Standard Model 
calculation is accepted, one finds [30]: 

a^(theory) = (11659180.5 ± 5.6) x 10"^°, (5) 
(experiment) = (11659208.0 ± 6.3) x 10"^°, (6) 

yielding a discrepancy [29] 

Aa^ = (27.5 ±8.4) x 10"^°, (7) 

which would be a 3.3-ct effect. In the plots discussed later, we display the corresponding 2-a 
range, namely 

10.7 X 10"^° < Attf, < 44.3 X 10"^°. (8) 



3.3 B Decay Observables 

We consider two constraints provided by B decay: one is the agreement between experiment 
and theory for b —>■ sj [31], and the other is the experimental upper limit on Bg —>■ yU"*"/^" 
decay. The recent measurements of t^u decay do not yet impinge significantly on the 

parameter space we explore in this paper. 

In the case of 6 — > 57, we use the estimate BR{b — > 57) = (3. 15 ±0.23) x 10~^ [32] for the 
SM contribution at NNLO 0, and the code of Gambino and Ganis E to calculate the MSSM 
contribution to the decay amplitude at NLO in QCD. As for the present experimental rate 
for 6 — > 57 decay, we use the range 

BR{b 57) = (3.55 ± 0.24+°;}^) x 10"^ (9) 

as recommended by the HFAG [33,34]. The first of the errors in (j9]) is the combined sta- 
tistical and systematic experimental error. The second set of errors result from theoretical 
uncertainties and corrections. These are combined linearly with the scale uncertainty in the 
calculation. We recall that b s'j joins in disfavouring /i < 0. 

In the case of Bg — > decay, we calculate the rate in the MSSM using [35,36], and 

we use the experimental upper limit 

BR{Bs /i+ fi-) < 1.0 X 10"^ (10) 

reported by CDF [37]. We also display in Figures H] - [5] projected future sensitivities of 
the Tevatron and LHC experiments (a factor of 5 times lower than the current limit). As 
already noted, the the impact of the Bg /U^ /i~ constraint is important at large tan /?, and 
increases as tua decreases and hence as Mj„ decreases. 

"^We note that the dominant theoretical error due to the renormahzation-scale uncertainty it is not 
Gaussian, and hence we add it hnearly rather than in quadrature with the other errors. 
^We thank Geri Ganis for a recent update to this code. 
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3.4 Neutralino Relic Density 



As aleady mentioned, we assume that the neutrahno LSP constitutes essentially all the cold 
dark matter, for which we consider the allowed range to be [20]: 

0.0855 < n^h^ < 0.1189, (11) 

as mandated by WMAP and other observations. 

As discussed in more detail in the Appendix, we have included in our calculation of 
the neutralino relic density, for the first time, all the processes for coannihilation between 
the three lightest neutralino states Xi,2,3, as well as with the lighter chargino and with 
sleptons. The importance of Xi~X2 — coannihilation has long been recognized within the 
context of the GUT-scale CMSSM [8,38]. Near the top-left boundary of the allowed region in 
the (/«i/2,''^o) plane, the lightest neutralino is Higgsino-like, near the bottom of the allowed 
region the lightest neutralino is bino-like, and the bino and Higgsino masses cross over along 
some intermediate contour. Near this cross-over line, and particularly where it intersects the 
left boundary of the allowed region in the (wii/2,'^o) plane, Xi ^ X2 ^ coannihilation is 
important in the GUT-scale CMSSM. 

In the GUT-less CMSSM, as we show later, there are interesting regions of the (mi/2, mo) 
plane at small Mj„ where the xs mass comes within (9(200) GeV of the X2 niass, and 
coannihilation processes involving the xs can no longer be neglected. The reason for this, 
despite the relatively large xs ~ X2 niass difference, is that the couplings of the Higgsino- 
like xs to relevant final states are significantly larger than the corresponding X2 couplings. 
Regions of the plane where X2 and xs are degenerate are present at most values of Mj„, 
though they typically occur when xi is much lighter than the other neutalinos. For low Mj„, 
however, there is in fact a near-degeneracy of all three of the lightest neutralinos as well 
as the lighter chargino. It is therefore necessary to include all coannihilations involving the 
three lightest neutralinos and the lighter chargino, as detailed in the Appendix. 

In addition, we implement here various improvements to our previous treatment of the 
dark-matter density in regions where rapid annihilation via a direct-channel Higgs pole is 
important. Specifically, we have included further crossed- channel contributions to WW, ZZ 
and less important processes. 

4 Evolving Impact of the Cold Dark Matter Constraint 

We now discuss the evolution of the dark matter constraint as the scale at which the soft 
supersymmetry breaking parameters are universal is lowered from the GUT scale. We assume 
rrit = 171.4 GeV in this analysis. Deviation by a few GeV from this value would result in 
some change to the exact positions and shapes of the regions preferred by WMAP, but our 
results are quite general. We recall that, as usual in the CMSSM, the value of the Higgs 
mixing parameter n is fixed by the electroweak vacuum conditions, leaving its sign as a free 
parameter. Motivated by and b —>■ sj, we consider only fi > 0, though a similar analysis 
could be carried out for negative /i. In Sections 14.11 and 14.21 we discuss in detail the effects 
of lowering the universality scale for two values of the ratio of the Higgs vevs, tan /3 = 10 
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and tan/? = 50. We take Aq = throughout Sections 14. II and and examine the impact 
of deviation from this assumption in Section I4.3[ Related mirage-mediation models are 
discussed in Section I4.4[ 

4.1 Low tan/3 

The evolution of the WMAP-preferred region in the {mi/2, rriQ) plane as the universality scale 
is lowered has been discussed previously in [15] for tan/5 = 10 and Mj„ > 10^^'^ GeV. The 
WMAP-preferred regions found in this analysis, along with constraints from colliders, are 
shown in Figs. [2] and [3] for several values of Mj„. To begin, we look first at the usual GUT- 
scale CMSSM scenario, shown in panel (a) of Fig. [2l One can see the x ~ ^ coannihilation 
region bordering the excluded stau LSP region for 330 ^ mi/2 < 900 GeV. Values of mi/2 
below this range are excluded by the LEP Higgs constraint. Near mi/2 = 900 GeV, the 
coannihilation strip dips down into the region where the f is the LSP. The focus point 
appears as a very thin strip tracking the border of the region excluded by the electroweak 
symmetry breaking condition at itlq > 1500 GeV. The LEP chargino bound also follows this 
boundary. The rapid-annihilation funnel is not present at tan j3 = 10 for Mj„ at the GUT 
scale, but will appear as the universality scale is lowered and also at larger tan (3. 

As found in [15], there are already changes as the universality scale is lowered to Min = 
10^^ GeV, shown in Panel (b) of Figure [21 The allowed focus-point region starts to separate 
from the LEP chargino bound, moving to larger mi/2- Notice also that this strip does not join 
smoothly with the coannihilation strip, but instead is deflected due to rapid h annihilation 
near mi/2 ~ 150 GeV. The region where the relic density falls in the WMAP range is thereby 
pushed inside the LEP chargino mass bound. However, this behavior occurs at low values 
of mi/2 which are excluded by the LEP Higgs bound as well. 

For Min = 10^'^ GeV, shown in panel (c) of Fig. [2], we notice that, as foreseen in Section O 
the regions excluded by the electroweak vacuum conditions and because the stau would be 
the LSP are encroaching further into the plane as Mj„ decreases, and the LEP Higgs bound 
is moving to larger mi/2. We see in panel (c) of Fig. [2] that the allowed focus-point region 
also dips further down, away from the electroweak vacuum condition boundary, while the 
coannihilation strip moves up and farther away from the region where the stau is the LSP. In 
fact, the focus-point and coannihilation regions connect, forming an slender atoll extending 
to {mi/2,mo) ~ (2850,2400) GeV (beyond the displayed region of the plane), inside which 
the relic density of neutralinos is too large. Another remarkable feature at this value of 
Min is the appearance of the rapid-annihilation funnel, familiar in the GUT-scale CMSSM 
at large tan/?, but an unfamiliar feature for tan/5 = 10. In the narrow space between the 
underside of the atoll and the thin WMAP-preferred strip lying approximately 100 — 200 
GeV below it, 2m^ ~ rriA and direct-channel annihilation processes are enhanced, causing 
the relic density to drop below the value determined by WMAP. 

As the universality scale is further decreased to Mj„ = 10^^'^ GeV, as shown in panel 
(d) of Fig. [2], the atoll formed by the conjunction of what had been the focus-point and 
coannihilation strips has shrunk, so that it lies entirely within the range of («^i/2, mo) shown 
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Figure 2: Examples of {mi/2, ttiq) planes with tan/? = 10 and Aq = but with different values 
ofMin. (a) The CMSSM case with Mi^ = Mqut ~ 2 x lO^^ GeV, (h) Mi^ = 10^^ GeV, (c) 
Min = 10^^ GeV and (d) = 10^^'^ GeV. In each panel, we show contours representing 
the LEP lower limits on the chargino mass (black dashed line), a Higgs mass of II4 GeV (red 
dashed), and the more exact (and relaxed) Higgs bound (red dot- dashed). We also show the 
region ruled out because the LSP would be charged (dark red shading), and that excluded by 
the electroweak vacuum condition (dark pink shading). The region favoured by the WMAP 
range VtcoMh'^ = 0.1045]'^o;qo95 has light turquoise shading, and the region suggested by — 2 
at 2-a has medium (pink) shading, with the 1-a contours shown as black dashed lines. 
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in panel (djj. We now see clearly two distinct regions of the plane excluded due to an excess 
relic density of neutralinos; the area enclosed by the atoll and the slice between the lower 
funnel wall and the boundary of the already-excluded f-LSP region. 

The four panels of Figure [3] show the consequences of lowering the universality scale even 
further, down as far as Mj„ = 10^ GeV. In panel (a) for Mj„ = 10^^ GeV, the focus-point 
and coannihilation regions are fully combined and the atoll has mostly filled in to become a 
small island of acceptable relic density. To the right of this island is a strip that is provided 
by the lower funnel wall. The strip curves slightly as mi/2 increases then takes a sharp 
plunge back down towards the boundary of the region where the stau is the LSP, a feature 
associated with the XX h + A threshold. Reduction in the universality scale from this 
point results in the lower funnel wall being pushed down into the excluded f LSP region and 
total evaporation of the island. 

As the universality scale decreases further in panels (b), (c) and (d) for Mj„ = 10^^ GeV, 
10^° GeV and 10^ GeV, respectively, we see only a small residual turquoise region at large 
mi/2 where the relic density is within the WMAP limits. At all other points in the visible 
part of the (^1/2, ""^o) plane the relic density of neutralinos is too low to provide fully the 
cold dark matter density preferred by WMAP et al. Of course, these SUSY models would 
not be excluded if there is another source of cold dark matter in the universe. 

In these last four panels, we notice that the boundary of the region where the stau is 
the LSP is retreating back down to smaller mo, as expected from the discussion of evolution 
with Min of the masses of the stau and the lightest neutralino given in Section O 

4.2 High tan/? 

The situation at larger tan/5 looks somewhat different at first glance. In the GUT-scale 
CMSSM case, shown in panel (a) of Fig. HI we see the familiar regions excluded because of a f 
LSP and the electroweak vacuum condition. The LEP Higgs and chargino mass bounds have 
impacts similar to those in the low-tan 13 scenario. The region excluded by 6 ^ 57 decay has 
grown substantially, and a new region excluded by the limit BR{Bs — > > 1 x 10~^ 

appears at low (mi/2,mo), which is, however, already excluded by other constraints. As 
for the relic density, the focus-point region is visible as a strip tracking the electroweak 
vacuum condition for mo > 1050 GeV, whereas the region preferred by WMAP is excluded 
by the LEP Higgs constraint at smaller mo. Along the excluded f LSP boundary, we see 
that the familiar coannihilation strip is truncated at low mi/2 by the Higgs and chargino 
mass constraints, and also by Bg fi~^ fi~ . Following this strip to larger mi/2, there is 
the familiar rapid-annihilation funnel, where 2m^j ~ m^ and the relic density is kept in 
the range preferred by WMAP by annihilations through the direct-channel A and H poles, 
which lifts away from the excluded region. 

However, at large tan/5, even small changes in the universality scale make a dramatic 
difference in the appearance of the regions preferred by WMAP. At Mj„ = 10^^'^ GeV, as 

^We note a string of bubbles intruding into the atoll, which are due to a significant enhancement of 
t— channel exchange in X2X2 h+ {H, A). The analysis of these possible regions of small relic density would 
require a complete treatment of poles, including finite-width effects, which we do not attempt here. 
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tan p = 10 , [I > , Mi„ = IQl^ GeV tan p = 10 , n > , Mj^ = lO" GeV 




mi/2 (GeV) mi/2 (GeV) 



Figure 3: Further examples of (mi/2,'"^o) planes with tan/5 = 10 and Aq = but with 
different values of Min: (a) Min = lO^^ GeV, (h) Mi„ = 10" GeV, (c) M,„ = 10^° GeV and 
(d) Min = 10^ GeV. The various contours and shadings are the same as for Fig. 2. 
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tan p = 50 , n > tan p = 50 , n > , Mi„ = lOl^-S GeV 




100 1000 2000 100 1000 2000 

mi/2 (GeV) mi/2 (GeV) 

Figure 4: Examples of {mi/2, mo) planes with tan/3 = 50 and Aq = but with different values 
of Mir,, (a) The CMSSM case with Mir, = Mqut ~ 2 x lO^^ GeV, (h) Mi^ = 10^^-^ GeV, 
(c) Min = 10^^ GeV and (d) Min = 10^^'^ GeV. In addition to the constraints enumerated 
in the caption to Fig. [1[ we also show the regions ruled out by b ^ s'y decay [31, 33, 34] 
(medium green shading) and black dot-dashed contours representing the current CDF limit 
on the rate of Bg —>■ fi~^ (1 x 10^^^ and a projected sensitivity of the Tevatron and the LHC 
experiments (2 x 10'^). 
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seen in panel (b) of Fig. HJ the coannihilation strip and rapid-annihilation funnel have joined 
to create a large funnel region that extends to (mi/2,mo) ~ (1850,2000) GeV. Inside the 
funnel boundary, effects such as rapid annihilation near the A pole and the coannihilations 
of neutralinos with light sleptons combine to cause the relic density to fall below the range 
preferred by WMAP. In this region of low relic density, the lightest neutralino is bino-like, 
and the dominant annihilations are into bb and rf pairs. As in the GUT-scale universality 
scenario, the focus-point region is cut off at mo ~ 950 GeV by the LEP Higgs constraint. At 
this universality scale, values of mi/2 > 600 GeV are compatible also with itlq > 2000 GeV, 
beyond the displayed region of the (mi/2,mo) plane. 

As the universality scale is further reduced to Mj„ = 10^^ GeV, we see in panel (c) of 
Fig. mthat the funnel is elongated further and opens wider at the top, while simultaneously 
the focus-point region falls significantly below the zone excluded by the electroweak vacuum 
conditions. In addition, the bulk region, where the upper funnel wall connects to the focus 
point, has now shifted to larger mi/2, so that it lies mostly outside the LEP Higgs bound. 
As in the other panels of this figure, the regions currently excluded by Bg fi~^ ix~ are also 
excluded by 6 ^ 57. We note that, as Mj„ decreases, the bulk and focus-point regions are 
moving to larger mi/2 more rapidly than the LEP Higgs constraint, resulting in a larger 
WMAP-preferred region at small mi/2 and rriQ. At the same time, however, the upper funnel 
wall is moving to smaller mi/2, causing the region between the focus point and the upper 
funnel wall (where the relic density is too large) to shrink. 

To illustrate how the relic density changes with mi/2 and its sensitivity to various inter- 
actions, we follow the evolution of the relic density for Mj„ = 10^^ GeV at a fixed value of 
mo = 1000 GeV. At very low mi/2 < 240 GeV, the electroweak symmetry breaking conditions 
would impose an unphysical solution for the weak scale value of the Higgs mass parameter, 
so this region of the plane is excluded, as discussed in Section [2J Near the boundary of the 
excluded region, /x < mi/2, so the LSP has a strong higgsino component and annihilations 
to light fermions keep the relic density low. As one moves to larger mi/2, the bino com- 
ponent increases, causing the relic density to increase accordingly, though it remains below 
the WMAP-preferred range. At mi/2 = 244 GeV, the XX W~^W~ threshold is reached 
and the relic density decreases dramatically, only to start rising again once the threshold is 
passed. By mi/2 = 280 GeV, the LSP has become bino-like, though it still has substantial 
higgsino components. 

Near mi/2 = 325 GeV, the relic density has risen to the range prefered by the WMAP 
measurements, and continues to increase until it exceeds the WMAP range. The thinness of 
the WMAP strip indicates the rate at which the relic density is increasing, reaching its peak 
value near mi/2 = 500 GeV. As mi/2 increases further, one approaches the broad {H,A) 
pole region, where s-channel annihilations cause the relic density to decrease dramatically. 
Thus, the upper funnel wall appears near mi/2 = 750 GeV, and the rehc density then 
continues to plummet until becomes large enough that the pole has been passed, at 
which point the relic density again increases until it falls within the WMAP range for a 
third time near mi/2 = 1080 GeV, forming the lower wall of the funnel region. As mi/2 
increases further, the relic density of neutralinos becomes too large to be compatible with 
the WMAP measurement. Near the border of the f LSP region, the relic density decreases 
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due to enhanced x — r coannihilations, however the effect is not sufficient to bring it down 
to the WMAP range. All values of mi/2 to the right of the lower funnel wall are excluded 
by the large relic density of neutralinos. 

When Min = 10^^'^ GeV, the focus-point region and upper funnel wall merge fully to 
form an island of acceptable relic density, extending from {mi/2,mo) ~ (400,850) GeV to 
large mo, parallel to the lower funnel wall, and with a width of ~ 200 GeV at its broadest 
point. 

In Fig. [5l as in the tan /3 = 10 scenario, we see the electroweak vacuum condition creep 
further down into the plane, as Min is further reduced. The f LSP region also retreats to 
smaller mo, because of the Min dependences of the sparticle masses discussed in Section El 
When the universality scale is Min = 10^^ GeV, as seen in panel (a) of Fig. [5l this island 
has submerged and disappeared as enhanced annihilations to bb and rf dominate even for 
2m^^ < TJiA- Coannihilations of Xo with Xi, where Xi denote the second- and third-lightest 
neutralinos, also play a significant role in the smallness of the relic density in this region. 
The only values of mi/2 and mo for which the relic density of neutralinos is in agreement 
with the WMAP measurement are in the thin strip that had been the lower funnel wall, and 
a narrow coannihilation strip adjacent to the f LSP region. To the left of the residual funnel 
strip, the relic density is below the WMAP value, whereas this value is exceeded in the 'vee' 
between the funnel and coannihilation strips at large mi/2. At Mj„ = 10^^ GeV, all values 
of mi/2 > 1230 GeV are excluded for mo < 2000 GeV. 

In panel (b) for Mj„ = 10^^ GeV, what is left of the lower funnel wall is also beginning 
to curve down. This is the same general behavior we observed in the tan/5 = 10 case. As 
the universality scale is slightly reduced, to Mj„ = 10^^'^ GeV (not shown), this strip bends 
down into the f LSP region at (mi/2,mo) ~ (2000, 1450) GeV. For Mj„ = 10^^ GeV, as seen 
in panel (c) of Fig. there remains only a small ellipse where the relic density falls in the 
region preferred by WMAP. The rest of the plane not excluded by the electroweak vacuum 
condition or the charged LSP constraint has a relic density of neutralinos smaller than that 
required by WMAP Q. This last remaining WMAP island evaporates as the universality scale 
is decreased to 10^^ GeV, as seen in panel (d) of Fig. [5l at which point the entire plane is 
disfavoured, in the sense that some additional source of cold dark matter would be required. 

4.3 Non-Zero 

To this point, we have considered all trilinear soft supersymmetry-breaking parameters to be 
zero at the unification scale, Aq = 0. Here we limit ourselves to a brief discussion of 7^ 
as preparation for the discussion of mirage-mediation models in the next Section. 

If Aq > 0, the RGEs generate correspondingly larger trilinear couplings at the weak 
scale. In addition, since the large loop corrections to /i depend on the values of the trilinear 
couplings, there is also an increase in /i. We therefore expect, based on the discussion 
in Section [21 that the region excluded by the electroweak vacuum condition decreases with 
increasing Aq. Other striking differences in the (mi/2, '^o) plane are in the constraints on the 

^We stress again that such regions are not excluded, provided there is another source of cold dark matter 
in the Universe. 
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tan p = 50 , n > , Mi„ = lOl^ GeV tan p = 50 , n > , M,,, = lO" GeV 
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Figure 5: Examples of {mi/2, mo) planes with tan (3 = 50 and Aq = but with different values 
ofMin. (a) Min = 10^^ GeV, (h) Min = lO^^ GeV, (c) Min = lO^^ GeV and (d) M^n = 10" 
GeV. The various contours and shadings are the same as for Fig. 4- 
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Higgs mass and the 6 ^ 57 rate. While the LEP Higgs constraint is dramatically relaxed 
for larger Aq, the region excluded by 6 ^ S7 increases in size, becoming the dominant 
constraint for low mi/2- Furthermore, since the off-diagonal elements of the squark mass 
matrix contain terms proportional to the negative of the trilinear couplings, when Aq is 
large these off-diagonal contributions can become large enough to drive the lightest stop 
quark mass below the LEP bound. As a result, we see a new excluded region emerge at low 
mi/2 and mo, where the lighter stop has m-j- < 220 GeV [39]. 

For Aq negative, the changes to the constraints discussed above are quite predictable. 
In this case, the RGE's generate correspondingly smaller weak scale trilinear couplings, 
resulting in a universally smaller /i. The LSP is then more more Higgsino-like over the whole 
plane. The LEP Higgs bound is strengthened, and the 6 — 57 rate becomes an insignificant 
constraint. 

The regions of the plane where the relic density of neutralinos is in the measured range 
also change shape for Aq 7^ 0. In general, these changes can be ascribed to one of two 
effects. First, in addition to the fx coannihilation strip, there may be an additional tx 
coannihilation strip, where the lighter stop is degenerate with the neutralino LSP. This 
feature is common in scenarios with large Aq and both the tx coannihilation strip and the 
excluded light stop region move further into the plane as Aq is increased. Secondly, we recall 
that the composition of the LSP depends on the ratio of n to Mi, the LSP being bino-like 
when Ml is small compared to /i and Higgsino-like if /i is small compared to Mi , as shown in 
panel (d) of Figure [H Since /z is enhanced everywhere in the plane when Aq > 0, we expect 
the LSP to be generically more bino-like than when Aq = 0. Similarly, we expect the LSP to 
be generically more Higgsino-like when Aq < 0. For ~ Mqut (not pictured), the LSP 
is strongly bino-like over most of the plane, so the main effects of 7^ are the above- 
mentioned modifications in the constraints, and the appearance of the tx coannihilation strip 
for large positive Aq. 

For lower unification scales, however, the LSP has more substantial Higgsino components, 
becoming Higgsino-dominated over much of the plane for very low Mj„. Larger fi means 
that the LSP will remain bino-like even for larger values of Mi, so in scenarios with Aq > 
the LSP is more bino-like and the heavier neutralinos with large Higgsino components are 
even heavier than when Aq = 0. These differences are clear at low Mj„, when the LSP is 
becoming Higgsino-like over much of the plane when Aq = but is still bino-like when Aq 
has a sufficiently large positive value. In panel (a) of Fig. [6l we show the (^1/2,^^0) plane 
for tan/? = 10, Mi„ = 10^^ GeV, and Aq = 1000 GeV. We note the similarity to panel (d) of 
Fig. [21 where Mj„ = 10^^'^. When Aq > 0, smaller values of /i appear only at values of Mj„ 
that are lower than in the Aq = cases previously discussed. In the same way, the Aq < 
case tends to mimic the effect of larger Mj„. With respect to the relic density of neutralinos, 
there is some degeneracy in the parameters Mj„ and Aq for regions of the (mi/2,mo) plane 
far from the tx coannihilation strip. 

We note that Aq oc M, where M = mi/2 or itlq, is also a viable possibility, the conse- 
quences of which, in light of the above discussion, are easily understood. In these cases, for 
small M, the plane will be similar to the Aq = case, while at larger M, the changes de- 
scribed above will be increasingly evident. A complete discussion of 7^ or non-universal 
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Aq is beyond the scope of this study. 
4.4 Mirage-Mediation Models 

Models in which supersymmetry breaking occurs through some combination of modulus and 
anomaly mediation are among those characterized by the apparent unification of gaugino 
and scalar mass parameters at an intermediate scale. As a result, these models have been 
termed mirage-mediation models [17] and the unification scale, the mirage messenger scale, 
is estimated to be ~ 10^° — 10^^ GeV. One distinctive feature of these scenarios is that the 
gaugino and scalar masses run both above and below the unification scale. Here, we discuss 
briefly the effect on our results of the additional running of the masses above the unification 
scale. 

The use of the RGEs to run the masses down from the input scale to the weak scale 
is unchanged, and the procedure for calculating the weak-scale observables is unchanged, 
regardless whether the soft supersymmetry-breaking mass parameters run above the unifi- 
cation scale. The chief difference derives from the fact that the value of fi is fixed by the 
electroweak vacuum conditions, which include a large dependence on the trilinear couplings 
as discussed in the previous section. When the trilinear couplings run from the GUT scale, 
becoming larger as the energy scale decreases, they attain larger weak-scale values than 
would be possible with running only below Min- Therefore, in mirage-mediation models /z 
receives a large contribution from the exceptionally large values of the trilinear couplings 
at the weak scale. The resulting picture for mirage-mediation models is similar to what 
one would expect from the GUT-less cases with Aq ^ 0, as discussed above. It should be 
noted that the trilinear couplings in mirage-mediation scenarios, as well as the other soft 
SUSY-breaking parameters, are specified at the GUT scale based on the particular mixture 
of modulus and anomaly mediation. The soft SUSY-breaking parameters are taken to be 
proportional to each other, with constants of proportionality determined by the modular 
weights and other considerations [17]. For simplicity, we consider only Aq = at the GUT 
scale. 

In panel (b) of Fig.[6]we show the (^1/2, "^o) plane with running of the gaugino and scalar 
masses both above and below the unification scale for Min = 10^^ GeV and tan /3 = 10. There 
is a broad region of acceptable relic density lying just above the excluded f LSP region. For 
comparison, in the standard GUT-less case for = 10^^ GeV shown in panel (b) of Fig. [3l 
as discussed already in Section 14. H the relic density of neutralinos is below the WMAP 2-a 
range throughout the plane, except in the small island just barely in view at mi/2 = 2000 
GeV. 

There are a few important differences worthy of note. First, the value of fi all over 
the plane is universally larger in the mirage-mediation scenario than in the cases discussed 
previously in this paper, which is attributed to the running of Aq from the GUT scale rather 
than Min- As a result, we expect the boundary of the region excluded by the electroweak 
vacuum conditions to be pushed back up into the upper left corner of the plane, as is seen. A 
second important consequence concerns the composition of the LSP. Recalling that the LSP 

^Such models are motivated, e.g., by the KKLT framework [19]. 
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Figure 6: Panel (a) shows the (mi/2,mo) plane for the GUT-less case with tan/3 = 10, 
Min = 10^^ GeV, and Aq = 1000 GeV. Panel (b) displays a scenario similar to that found in 
mirage-mediation models, where the soft supersymmetry- breaking parameters are universal 
at Min = lO^"*^ GeV, but run both above and below this scale. The weak-scale values of the 
neutralino and chargino masses, as well as the pseudoscalar Higgs mass and ^, are shown in 
panel (c) for the usual GUT-less case with = and Min = ^0^^ as shown in panel (b) of 
Fig. Panel (d) shows the same information as panel (c) for the mirage-mediation case. 
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is bino-like as long as Mi is much smaller than /i. The fact that n is larger in the mirage- 
mediation case implies that the cross-over when Mi ^ /i takes place at a lower unification 
scale than was found in panel (d) of Fig. [H In fact, the LSP is bino-like over most of the 
plane in the mirage-mediation case shown in panel (b) of Fig. [6l whereas it has large Higgsino 
components for much of the standard GUT-less plane for the same value of Mm. Similarly, 
the heavier neutralinos, which have large Higgsino components, are even heavier due to the 
enhancement in /i in mirage-mediation models. This effect can be seen clearly by comparing 
panels (c) and (d) of Fig. [61 

5 Neutralino-Nucleon Cross Sections 

Direct searches for dark matter particles such as the Cryogenic Dark Matter Search (CDMS) 
[40] and other experiments look for evidence of weakly-interacting massive particles (WIMPS) 
through scattering on nuclei. In this section, we present the predictions for neutralino- 
nucleon scattering cross sections in the scenarios discussed above [41] - [48]. 

The low-energy effective interaction Lagrangian for elastic x-iiucleon scattering can be 
written as 

L = a2tXl^i^XQtlfJil^(lt + asiXXQiQi , (12) 

where terms that make velocity-dependent contributions to the cross section have been ne- 
glected, and the constants a2i and a^i are defined as in Ref. [43]. In computing the scalar 
cross section, we have assumed the pi-nucleon E term to be 64 MeV (see [49] for the sensi- 
tivity of the elastic cross section to this assumption). Summation over the quark generations 
is implied, with up- and down-type quarks labeled by the subscript i. The cross section 
can be broken into a spin-dependent part arising from the term proportional to a2i and 
a spin-independent (scalar) part from the term proportional to a^i. The spin-dependent 
cross section is, in general, larger than the scalar cross section. However, since the whole 
nucleus participates coherently in spin-independent interactions, it is primarily the scalar 
cross section that is probed by current direct-detection experiments. On the other hand, the 
spin-dependent scattering cross section on the proton plays an important role in the capture 
and annihilation rates inside the Sun. 

Figs. [7] and [8] show scatter plots of the spin- dependent and scalar cross sections for elastic 
X-nucleon scattering. We plot the cross sections as functions of the neutralino mass for points 
in the {mi/2, mo) plane where the relic density of neutralinos is less than the 2-a upper limit 
from WMAP (as first examined in Ref. [50]) with the assumption of universality at the GUT 
scale relaxed. For the cases where the relic density is smaller than the central WMAP value, 
indicating that there must be another source of astrophysical cold dark matter, we plot the 
cross section scaled by the ratio of the relic density of neutralinos to the central density of 
cold dark matter inferred from WMAP measurements of the CMB. These results can be 
compared with the direct-detection limits available from CDMS and other experiments. In 
each figure, we also show the CDMS II limit for the scalar part of the neutralino-nucleon 
cross section [51]. Current limits on the spin-dependent cross section are cr^„ < 10^^ pb [52], 
which lies outside the range we have plotted in Figures [7] and [HI We require that the 
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lightest neutralino be the LSP and that electroweak symmetry be broken, as usual. The 
LEP constraint on the chargino mass has been applied, as discussed in Section 13.11 Different 
colors in Figures [7] and M indicate whether the point lies within the region excluded by 
b ^ sj decay or the LEP Higgs mass constraint. The dark blue (striped) regions are the 
spin-dependent (scalar) cross sections that pass all these constraints. Lighter (green) regions 
in each panel fail the relaxed LEP Higgs constraint. At large tan/3, when the constraint on 
the rate of b —>■ sj becomes dominant, we show in red the regions that fail this constraint 
but pass all others. 

The results for tan/3 = 10 are shown in Figure [7] for four different values of Mj„. The 
spin-dependent cross sections lie above the scalar cross sections in each of the four panels. 
As Min is lowered, the number of points increases dramatically and they spread to larger 
m^. This is due to the fact that the rehc density over all of the (mi/2,'n^o) plane decreases 
as Min decreases, so that less and less of the plane is excluded by having an excess relic 
density [^. 

We turn our attention first to the usual GUT-scale CMSSM, in which the relic density 
is too large over most of the {itli/2,itlq) plane. Within the allowed regions for each of the 
spin-dependent and scalar cross sections, we can identify two separate behaviours. First, 
there is a region stretching out to ~ 350 GeV where the cross section may vary over 
as much as an order of magnitude for some values of m^. This feature corresponds to the 
coannihilation strip, which is shown in Figure [2] to dip into the f-LSP excluded region near 
^1/2 = 900 GeV. The variation in the cross section in this coannihilation strip region at low 
in panel (a) is due to the separation of the coannihilation strip from the boundary of 
the f-LSP region at low mi/2. The cross sections for points lying between the coannihilation 
strip and the forbidden f-LSP region, where the relic density of neutralinos is too low, are 
scaled down to reflect the fact that in these cases the neutralinos can provide only a small 
fraction of the cold dark matter in the Universe. 

The second region lies within 80 GeV < < 170 GeV. In the case of the spin-dependent 
cross section, the cross sections in this region are clearly separated from those due to the 
coannihilation strip. This second region of acceptable cross sections comes from the focus- 
point region which, for Mj„ = Mquti occurs at large mo and small mi/2- It should be noted 
that, if we were to consider values of uiq > 2000 GeV, the focus-point region would extend 
to larger (mi/2,mo), so analogous focus-point cross sections would extend also to larger m^. 
In the focus-point region, the fact that the lightest neutralino acquires substantial Higgsino 
components leads to an enhancement in the spin- dependent cross section due to Z exchange. 
Simultaneously, the scalar cross section becomes dominated by neutral Higgs exchange as 
the neutralino becomes Higgsino-like. 

Panels (b), (c), and (d) show the neutralino- nucleon cross sections as functions of the 
neutralino mass for Mj„ = 10^^, 10^^, and 10^° GeV, respectively. The changes in the cross 
sections as Mj„ is lowered may be understood by referring to the corresponding (mi/2,mo) 
planes from Figures [2] and [31 When = 10^^ GeV, the focus-point region becomes more 
prominent, separating from the boundary of the region excluded by the electroweak vacuum 

^In fact, for Min = 10^" GeV as shown in Panel (d), the constraint on the relic density does not exclude 
any points, but serves only as a scale factor for the cross sections. 



21 




(GeV) (GeV) 



Figure 7: Neutralino-nucleon cross sections as functions of the neutralino mass for tan (3 = 10 
andAo = but with different values of Min. (a) Min = Mgut ~ 2x10^^ GeV, (h) Min = 10^^ 
GeV, (c) = 10^2 GeV and (d) M,„ = lO^o GeV. 
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Figure 8: Neutralino-nucleon cross sections as a function of neutralino mass for tan (3 = 50 
and Ao = but with different values of Min. (a) Min = Mqut ~ 2x 10^^ GeV, (h) Min = 10^^ 
GeV, (c) = 10^^ GeV and (d) M,in = lO^^ GeV. 
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condition. In fact, for the portion of the (mi/2,'^o) plane shown in Fig. [21 the focus-point 
region extends to larger than the coannihilation strip. The two regions are seen as merged 
in panel (b-d). 

As we proceed to panel (c), most of the (mi/2,mo) plane results in a rehc density of 
neutralinos that is within or below the cosmologically-preferred range. As a result, there is a 
uniform distribution of possible cross sections up to ^ 650 GeV. The upper boundaries 
of the scalar and spin-dependent cross sections in panel (c) come from regions in the plane 
where the relic density is largest and mo is lowest, i.e., from the WMAP-preferred regions 
found at low mo. The continuous WMAP region that extends from the f-LSP boundary 
to larger mo and mi/2 is responsible for this uniform upper limit for the cross sections for 
m^ < 650 GeV. Near mi/2 = 1100 GeV, however, a new region of preferred relic density 
emerges at lower mo, leading to a bump in the neutralino-nucleon cross sections that extends 
to the largest values of m^ considered here. 

This same behavior is observed in panel (d), where Mj„ = 10^*^ GeV. The relic density of 
neutralinos falls within the WMAP range only in a small region of the (mi/2; ''^o) plane with 
mi/2 > 1700 GeV and is too small elsewhere, but similar increases and decreases in the relic 
density where different annihilation channels dominate are evident. We also point out that, 
since the relic density is lower than the WMAP range over most of the plane for Mj„ = 10^^ 
and 10^° GeV, we clearly see the maximum weak-scale value of m^, which corresponds to 
mi/2 = 2000 GeV, decrease between panel (c) and panel (d). For Min < 10^^ GeV, the LSP 
becomes Higgsino-like, with m^ ~ /i and fi decreasing rapidly as is lowered, as discussed 
in Section [21 For tan/5 = 10, the cross sections excluded by CDMS come only from points 
in the (mi/2,mo) plane that also fail the relaxed LEP Higgs constraint. 

In Figure [8] we show the neutralino-nucleon cross sections for tan (3 = 50 with Mj„ = 
Mgut, 10^^, 10^^ and 10^^ GeV. Although the cosmologically-preferred regions of the (mi/2, 
mo) plane are somewhat different from those for tan/? = 10, the plots in Fig. [H] look qual- 
itatively similar to those in Fig. [3 In panel (a) there is a clear separation between the 
cross sections from the focus-point region and those from the coannihilation strip and the 
beginning of the rapid-annihilation funnel. Since the funnel region of acceptable relic den- 
sity pictured in panel (a) of Fig. [H extends to mi/2 ~ 1850 GeV, we find values of the cross 
sections out to rriy ^ 850 GeV. We note that some of the scalar cross sections for m^ < 200 
GeV that pass all other constraints outlined above have been excluded by CDMS. 

In panel (b), where M^n = 10^^ GeV, the two regions are still distinct. The lower bulk 
of cross sections comes now from points inside the fully-developed rapid-annihilation funnel, 
seen in panel (c) of Fig. [H Again, we note that had we extended our analysis to larger values 
of mi/2 and mo, acceptable cross sections would be found also at larger m^. 

When Min = 10^^ GeV, shown in panel (c), the upper funnel wall has passed through 
the focus point, and only the lower funnel wall remains. Regions to the left of this wall in 
the (mi/2,''^o) plane are essentially inside the funnel and have a very low relic density of 
neutralinos, whereas the relic density is too large to the right of the wall. Consequently, we 
see in panel (c) that, at low m^, the scalar and spin-dependent cross sections span several 
orders of magnitude. 

As in the case when tan/5 = 10, at low the relic density of neutralinos falls within 
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or below the WMAP range over all of the (mi/2, mo) plane, so none of the plane is excluded 
by the constraint on the relic density. This is the case in Panel (d) , where Mj„ = 10^^ GeV. 
The situation remains unchanged as Mj„ is further decreased. 

We note that the scalar cross sections are generally larger at large tan/3. In fact, some 
of these cross sections are already excluded by ZEPLIN-II as well as CDMS II, which both 
probe WIMP-nucleon scalar cross sections as low as a few xlO~^ pb [51,53]. A sensitivity 
of 10~^ pb for ^ 100 GeV is expected for SuperCDMS Phase A with seven towers 
deployed [54]. Many direct dark-matter search experiments plan to use Xenon or Argon 
as an alternative target material for which the sensitivity scales linearly with the detector 
mass. The Argon Dark Matter exeriment (ArDM) expects to probe spin-independent cross 
sections as low as 10~^° pb with a one-tonne detector operating for one year [55]. Results 
from direct detection experiments will provide a useful complement to searches for SUSY 
signatures at colliders. 

6 Summary 

We have examined the impact of lowering the scale of unification of the soft supersymmetry- 
breaking parameters of the CMSSM on phenomenological, coUider and cosmological con- 
straints. In order to carry out this study, we accounted for coannihilations involving the 
three lightest neutralinos, the lighter chargino, and relevant sleptons and squarks. We ex- 
plored tan (3 = 10 and tan (3 = 50, Aq ^ 0, and a specific case similar to those found 
in mirage-mediation models. Intermediate unification scales result in the appearance of a 
rapid-annihilation funnel even at low tan (3, and the merging of this funnel and the focus- 
point region as Mj„ decreases. As the unification scale is lowered below a critical value 
dependent on tan/5 and other factors, the relic density of neutralinos becomes too low to 
account fully for the required relic density of cold dark matter over all or nearly all of the 
(mi/2, m-o) plane. These values of Mj„ are disfavored in the sense that there must be another 
source of astrophysical cold dark matter in the universe. 

We have also presented the neutralino-nucleon cross sections for several values of Mj„ 
at tan f3 = 10 and tan j3 = 50. We find that the spin-independent neutralino-nucleon cross 
sections for regions of parameter space favored by cosmology are beginning to be exchidcd 
by CDMS and other direct detection WIMP searches, although viable cross sections span 
several orders of magnitude. We look forward to stronger limits on the spin-independent 
cross sections as direct-detection WIMP searches become more sensitive in the near future. 

The analysis in this paper has shown that lowering the scale of unification in even the 
simplest CMSSM model may alter significantly the phenomenological expectations for both 
collider and non-collider experiments. It has also revealed novel effects in the calculation of 
the relic neutralino density, such as the importance of multi-channel neutralino and chargino 
coannihilation processes. However, we have done little more than scratch the surface of 
possibilities since, for example, we have not considered in detail scenarios with different 
values of Aq, let alone non-CMSSM scenarios or more realistic mirage-mediation models. 
Another interesting and important question for the future is the accuracy with which the 
effective unification scale could be estimated on the basis of future collider experiments. We 
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hope that this work will trigger future studies of these and other related issues. 
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A Neutralino and Chargino Coannihilations 

In most standard CMSSM scenarios, the LSP is a bino-like neutralino in many of the regions 
of parameter space relevant to cosmology, possibly with a significant Higgsino admixture. 
When the relic density falls near the range favoured by WMAP and other measurements, 
the neutralinos are typically not degenerate, and therefore there is no opportunity for coan- 
nihilations of the LSP with other neutralinos, or with charginos, to bring the relic density 
down into the range preferred by cosmology. The only case in which it is necessary to in- 
clude coannihilations involving neutralinos and charginos occurs when the neutralino LSP is 
Higgsino-like, which arises when /i < Mi^a situation that may arise at large mo in the focus- 
point region of the GUT-scale CMSSM H. In such it is possible for the lightest and 
second-lightest neutralinos to be degenerate with each other and with the lightest chargino. 
Thus, at large tuq and small mi/2 in the GUT-scale CMSSM, coannihilations between the 
lightest and second lightest neutralinos and with charginos must be included 0. 

In GUT-less scenarios, the lightest neutralino becomes Higgsino-like at low Mj„, as dis- 
cussed in Section |2], so it is necessary to include coannihilations involving the two lightest 
neutralinos and the lightest chargino as discussed above. However, there is also a region 
of parameter space where additional coannihilations become significant. When the LSP is 
mixed and nearly degenerate with the second lightest neutralino and the lightest chargino, 
in some circumstances the third-lightest neutralino may also be nearly degenerate. 

In panels (a) and (b) of Figure [9] we show the masses of all neutralinos and charginos 
as functions of Mj„ for two different points in the plane. We recall from the 

discussion in Section [2] that the slope of the curve describing the LSP mass is an indication 
of its composition: when the neutralino mass increases as Mj„ decreases, it is gaugino- 
like, and when it decreases as Mj„ decreases, tracking it is Higgsino-like. Panels (a) 
and (b) show that the masses of the LSP, the second lightest neutralino and the chargino 
are nearly degenerate when the LSP is Higgsino-like, indicating the necessity of including 
coannihilations involving all three states. Moreover, just at the point where the LSP changes 
from bino-like to Higgsino-like, the mass of the third-lightest neutralino dips down near the 
masses of the two lighter neutralinos and the chargino. Panels (c) and (d) compare the 
masses of the neutralinos and charginos as functions of mo for the CMSSM case with GUT- 
scale universality and Mj„ = 3 x 10^^ GeV for fixed mi/2 = 1000 GeV. One can see in panel 
(c) that, for the GUT-scale case, there is no degeneracy of the LSP with other neutralinos 

^"^This situation may also occur in some models with non-universal Higgs masses. 

^^It should be noted that f-x coannihilations are known to be of general importance in the CMSSM, since 
they give rise to the so-called coannihilation strip. 
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or charginos. The LSP is strongly bino-like, and therefore its mass is related to mi/2, as in 
([2]), with only a very weak dependence on rriQ through higher-order corrections. The same 
scenario is shown in panel (d) for Mj„ = 3 x 10^^ GeV. In this case, however, we see that 
mass degeneracies are apparent over a wide range of values of rriQ. 

For some values of Mj„, the near-degeneracy of the third-lightest neutralino with lighter 
neutralinos and the chargino occurs precisely where the relic density of neutralinos is near 
the cosmologically preferred value. For example, for Mj„ = 10^^ GeV and tan/3 = 10, the 
shape and location of the WMAP strip running through (1500, 1000) GeV can shift by as 
much as 200 GeV in mo if coannihilations are not properly included. 

Thus, the calculations of the coannihilation processes that previously were included for 
the second-lightest neutralino have here been calculated also for the third-lightest neutralino, 
including those of X2 with ^3. Table [1] shows the initial states for all the calculated annihila- 
tions and coannihilations of neutralinos and charginos used in the analysis here. In addition 
to those outlined below, coannihilations of all these neutralino and chargino species with 
sfermions were calculated, as well as the corresponding sfermion-antisfermion annihilation 
processes. 



XiXi 


X1X2 


XiX+ 


X2X2 


X2X3 


X2X+ 


X3X3 


XsXi 


X1.X+ 




X+X- 





Table 1: Initial states of interactions included here in the calculation of the relic cold dark 
matter density, where xi is the LSP and X2(3) is the second (third) -lightest neutralino. 

As noted in the main text, we take the opportunity in this paper to improve on our 
previous treatment of the rapid- annihilation region and to correct certain coding inaccuracies 
which, however, have no visible effects on the results we present. 
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